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ABSTRACT 

We report on an outburst of the high mass X-ray binary 4U 0115+63 with a pulse period of 3.6 s in 2008 March/April as observed 
with RXTE and INTEGRAL. During the outburst the neutron star's luminosity varied by a factor of 10 in the 3-50 keV band. In 
agreement with earlier work we find evidence for five cyclotron resonance scattering features at ~10.7, 21.8, 35.5, 46.7, and 59.7 keV. 
Previous work had found an anti-correlation between the fundamental cyclotron line energy and the X-ray flux. We show that this 
apparent anti-correlation is probably due to the unphysical interplay of parameters of the cyclotron line with the continuum models 
used previously, e.g., the negative and positive exponent power law (NPEX). For this model, we show that cyclotron line modeling 
erroneously leads to describing part of the exponential cutoff and the continuum variability, and not the cyclotron lines. When the 
X-ray continuum is modeled with a simple exponentially cutoff power law modified by a Gaussian emission feature around lOkeV, 
the correlation between the line energy and the flux vanishes and the line parameters remain virtually constant over the outburst. 
We therefore conclude that the previously reported anti-correlation is an artifact of the assumptions adopted in the modeling of the 
continuum. 
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1. Introduction 

In all accreting X-ray pulsars the accretion geometry is strongly 
affected by the magnetic field of the neutron star, which has 
values of the order of a few 10 12 G at the magnetic poles. 
The accreted matter is forced to follow the magnetic field lines 
from the Alfven radius inward to the accretion columns, lo- 
cated at the magnetic poles of the neutron star. The shape of 
these columns remains an open question. The simplest models 
assume cylinder or slab-like geometries, although more compli- 
cated shapes have been discussed (Meszaros 1984[|Becker et al.| 



2012). On each magnetic pole of the neutron star a hot spot 
is created, emitting thermal photons into the accretion column 
( |Davidson||1973| ). These photons are upscattered in energy by 
inverse Compton scattering within the accretion column, form- 
ing the hard X-ray spectrum ( |Becker & Wolff 2007 
et al.|2009| ). Depending on the source luminosity, and t 



Ferrigno 



iius mass 

accretion rate, M, and deceleration mechanism, a radiation dom- 



inated shock can form above the surface of the neutron star. See 
Bec ker et al.| ( |2012| ) for a discussion of the different modes of 
accretion in the accretion column. 

Due to the strong magnetic field of accreting X-ray pulsars, 
cyclotron resonance scattering features (cyclotron lines; CRSFs) 
are observable in many X-ray pulsar s. These features have been 
reported for about 20 X-ray pulsars (Caballero & Wilms 2011 ). 
CRSFs are produced by photons generated in the ac cretion col- 
umn of the neutron star interacting with electrons ( Schonherr 
|et al. 12007} |Meszaros| 1 992, and references therein). The motion 
of these electrons perpendicular to the 5-field is quantized into 
Landau levels with energy differences of about 
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Photons inelastically scattering off these electrons will imprint 
absorption-line like features onto the spectral continuum. 
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Fig. 1. RXTE/ ASM lightcurve of 4U 0115+634 (2-10 keV). The 
inset shows a close-up view of this lightcurve, showing the 2008 
outburst in detail. Arrows indicate the times of pointed RXTE 
and INTEGRAL observations. Numbers indicate the observation 
epochs used in the data analysis (see also Table [TJ. 



Measuring the energy of the CRSF yields information about 
the 5-field strength in the line forming region. Since to first or- 
der the 5-field of a neutron star can be assumed to be dipole- 
like, in many models the centroid energy of CRSFs depends on 
the characteristi c emission height, i.e., the altitude of the line 
forming region ([Basko & Sunyaev||1976| |Burnard et al.||199T| 
Mihara et al. 2004J). Since this altitude depends on the complex 
relation between the deceleration of the accreted material and 
the mass accretion rate, it is expected to depend on the source 
luminosity. Such a dependence has indeed been observed. Three 
types of behaviour have been observed. One type shows a pos- 
itive correlation between the X-ray luminosity Lx and the cy - 
clotron line energy E cyc (e.g., Her X-l, Staubert et al. 2007). 
The other type sh ows a negative correlat ion between L x and E cyc 
(e.g., V0332+53, |Nakajima et al.|2010| and references therein). 
Finally, there are also cyclotron sources known with a constant 
cyclotron line energy (e.g., A0535+ 26,|Caballero et al. 200 7]). 

As suggested by Sta ubert et al.|(|2007|) and Klochkov et al. 
( |2011| ) and then quantified by |Becker et aL] ( |2012| ), these differ- 
ent relations between L x and E cyc could be due to the different 
deceleration mechanisms for the accreted material. Above a cer- 
tain critical luminosity, L cr i t , the accreted material is decelerated 
mainly by the radiation pressure. In this situation, an increase in 
luminosity stops the material earlier, leading to the line forming 
region being in a region of lower 5-field strength and yielding 
a negative correlation between E cyc and Lx. For very small lu- 
minosities, the deceleration is dominated by the gas pressure. In 
this regime larger mass accretion rates will compress the accre- 
tion mound and thus lead to a positive correlation between the 
luminosity and the cyclotron line energy. Finally, at intermediate 
luminosities one expects a regime where the accretion column is 
fairly stable and the deceleration is dominated by Coulomb in- 
teractions. 

4U 0115+634 ( [Giacconi etal.|1972j | is one of the X-ray pul- 
sars for which CRSFs have been studied in great detail (see, 



Table 1. Log of observations for the 2008 outburst of 
4U 01 15+634. The first part of the table contains the RXTE ob- 
servations, the second part the ISGRI data. 
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Notes. For RXTE, the first column contains the number of the Obs- 
ID after "93032-". For ISGRI, it contains the revolution number. b For 
RXTE, the two numbers correspond to the PCA and HEXTE exposures, 
respectively. (c) Epoch for data grouping, see text for details. 
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The source shows 3.61 s pulsation and has an or- 
bital p eriod of -24.3 d (Comin sky et al.||1978[|Rappaport et aT 



1978 



et al. 



The optical counterpar t is the 09e star V635 Cas ( [Johns 
T978}|Unger et al.|1998|). 



In previous outb ursts, CRSFs h ave been detected up to the 
fifth harmonic (Heind i~et al.||1999] ). This is the largest number 
of detected CRSFs in an accreting X-ray pulsar. In later work, 
a negative correlation between Lx and E C y C was found (e.g ., 



Nakajim a etal.|2006| [Tsyganko v et al.|2007[ [Miiller et al.|2QTQ] ) 9 
meaning that the source should be located in the supercritical lu- 
minosity regime. Other authors even found a discontinuous be- 
havior of the fundamental line energy and explained this by a 
sudden change of the accretion conditions ( |Li et al.|2012| ). 

In this work we revisit the behavior of the CRSFs in 
4U 0115+634. Focusing on data taken during a giant outburst 
in 2008 spring (Fig. [TJ. We show that the previously claimed 
relation between the souce flux F x and E cyc strongly depends 
on the choice of the underlying continuum model. In Sect. [2] we 
summarize the observations. In Sect. Owe describe the differ- 
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ent continuum models and the flux and time dependence of the 
continuum and cyclotron line parameters. We argue that the cy- 
clotron line is only modelled correctly when using an exponen- 
tially cutoff power law continuum with a strong emission feature 
around lOkeV. Furthermore, alternative continua models, which 
have been used in some earlier investigations, result in incorrect 
cyclotron line parameters which can cause the apparent relation 
between source flux Fx and fundamental cyclotron line energy 
Eq. We summarize this paper in Sect. [5] 



2. Observations and data reduction 

The outburst of 4U 0115+634 started 2008 March 12, after a 
phase of almost four years of quiescence, when Swift/TSAI de- 



tected a significant increase of the X-ray flux (Kri mm et al 
2008). The outburst lasted about 40 days and exceeded a flux 



of ~200mCrab in the 2-12 keV RXTE/ ASM band (Fig.[l}. The 
whole outburst was monitored with pointe d observations by 
the Rossi X-ray Timing Explorer (RXTE, |Bradt et al.|[T993] ) 
and several pointed observations by the INTErnation al Gamma - 
Ray Astrophysics Laboratory (INTEGRAL, \Win]der et al.|2003) . 
Here, we use data from RXTE's Proportional Counter Array 
(PCA, (Jahoda et al.]|2006]) and the High Energ y X-Ray Timing 
Experiment (HEXTET lRothschild et aLl^998] ), as well as the 



INTE GRAL Soft Gamma-Ray Imager (ISGRI, |Lebrun et aL 
2003 ). Data were reduced with our standard analysis pipelines, 
based on HEASOFT (v. 6.11) and INTEGRAL OS A v. 9.0. 
Data modeling was performed with the newest version of 



the Interactiv e Spectral Interpretation System (ISIS, Houck & 
|Denicola|2000l ). 

The PCA consisted of five proportional counter units 
(PCUs), which were sensitive betweeen 2 and 60keV. Since 
PCU 2 is known to have been the best calibrated one ( |Jahoda| 
et al.|2006| ), we only extracted spectra from the top layer of this 



PCU in the standard2f mode. The hard X-ray spectrum from 
1 5-250 keV was monitored with HEXTE. HEXTE consisted of 
two independent clusters, A and B, both arrays of four phoswich 
scintillation detectors. The "rocking mechanism" of these clus- 
ters allowed a near realtime measurement of the background. 
Since this mechanism failed for cluster A before our observa- 
tions, we only used data from cluster B. The other instrument 
for the hard X-rays was the pixelized CdTe dector ISGRI be- 
hind a Coded Mask on board the INTEGRAL satellite. It covers 
the energy range between ~18 keV and 1 MeV. We extracted all 
available ISGRI spectra, for which 4U 0115+634 was less than 
10° off-axis. Table[T]contains an observation log (see also Fig. [I]). 

For the spectral analysis we used data between 3 and 50 keV 
for PCA, and between 20 and lOOkeV for HEXTE and for 
ISGRI. A systematic error has been added in quadrature to the 
PCA and ISGRI data, us ing the canonical values of 0.5% and 
2.0% ( [Jahoda et 81^20 06 , and IBIS Analysis User Manuaf}, re- 
spectively. To improve the signal to noise ratio of the spectra, 
we averaged the data over 14 data blocks in time. These data 
blocks are defined in the fifth column of Table Q] and shown in 
Fig. [T] We omitted one ISGRI dataset from our analysis (rev- 
olution number 677), because no simultaneous RXTE data are 
available at that time and the source was almost in the off state 
during that observation. 



3. Spectral analysis 

3.1. Pulsar Continuum Models 

A large problem when analyzing X-ray spectra from X-ray bi- 
nary pulsars is to find a good description of the broadband con- 
tinuum. In most cases, semi-physical models, which are some 
variant of a power law with a high energy exponential cutoff 
(Kreykenb ohm et aL]|2002| and references therein), are used to 
describe the spectra. These models can be justified by consider- 
ing the Comptonization of soft photons in an accretion column 
(Becker & Wolff 2007, and references therein). The most basic 
exponentially cutoff power law model, called Cut of fPL in ISIS 
and XSPEC flArnaud|1996] ), has the form 



Cutof fPL(E) oc E~ T exp (-E/E Md ) , 



(2) 



with the photon index , T, and the folding energy, £f old (see, e.g., 
Schonher r et al.|20 07). However, this simple model is often, but 
not always, insufficient to describe the continua of observed X- 
ray pulsars, as the curvature in the continuum is often seen to 
start at different energies. To obtain a good spectral fit, therefore, 
models with more free fit parameters and more complex shapes 
are sometimes required. The most simple of these models is the 
power law with high energy cutoff 



PLCUT(£) oc 



\E T for £ < £ break 

[E~ r exp ((£ b reak - E)/E Md ) for E > £ break 

(3) 

where for typical X-ray pulsars £ brea k ~ lOkeV and which has 
an unphysical "kink" at £t>reak ? which can result in line-like resid- 
uals around this energy ([Kretschmar et al. |1 997 [ [Kreykenbohm 



et al. 



et al. 



1999|). For this reason Klochk ov et al.| (2008 ) and Ferrigno 



(2009) smooth the continuum around this break, using a 



model of the form 

IE~ T for E < £ brea k - A£ 

AE 3 + BE 2 + CE + D for \E - £ break | < A£ 
E~ Y exp (-E/E Md ) for E > £ break + A£ 

(4) 

where the coefficients A, B, C, and D are chosen such that the 
continuum and its first derivative are continuous at E = £t>reak ± 
AE and where AE ± but sufficiently small. 

Other common choices to avoid the break are the so called 
Negative and Positive power-law Exponential (NPEX, see, e.g., 
Makishi ma et al.|[l999| and references therein), a continuum 
model with five free parameters of the form 

NPEX(£) oc (E~ Tl + aE +T2 ) exp (-E/E Md ) , (5) 

with photon indices Y\ 2 > and the "Fermi Dirac Cutoff" 
(FDCUT, see |Tanaka|1986| ) 



FDCUT(£) oc E~ T 



1 



1 + exp ((E - £ b reak)/£fold) 



(6) 



Common to all of these empirical continua is that despite the 
smoother transition around £ brea k they often need to be modified 
further by a broad excess or depression around lOkeV. This so- 
called "lOkeV feature" has been observed in many X-ray pulsar 
spectra, either i n absorption or in emission, although its origin is 
still not clear (Coburn 2001). This feature is usually modeled as 
an additive Gaussian line, with the centroid energy Eq, the width 
<tg, and the flux Aq. 



http://www.isdc.unige.ch/integral/analysis#Documentation 



Often one fixes T 2 = 2, see, e.g., Nakajima et al. (2006 1. 
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In addition to those empirical models, based on the theoret- 
ical description of bulk and thermal Comptonization in mag- 
netized accretion colums of Becker & Wolff ( 2007 , and refer- 
ences therein), Ferrigno et al.| ( |2009| ) have developed a contin- 
uum model depending on physical parameters like the culumn 
radius, the accretion rate, the electron temperature, the #-field 
strength. A hybrid model which combines the broadband con- 
tinuum with CRSFs is currently in preparation ( [Schwann et al.| 
2012 ). These physics-based continuum models are still in active 



development and, because of the large number of free parame 
ters, only suited for very high signal to noise data. For this reason 
most observational work uses one of the continua described by 

For sources with cyclotron lines, the continuum model is 
modified with a description of the cyclotron line. Since phys- 
ically self consistent line models such as those described by 
|Schonherr et al.| ( |2007| ) are still not available for all sources, em- 
pirical line models are generally used. The two most common 
CRSF models are those using a Gaussian or a pseudo-Lorentzian 
optical depth profile (see discussion by Nakaji ma et al.||20T() ). 
Here, we model the CRSFs with a pseudo-Lorentzian optical 
depth profile (Miha raet al.|[P9 90), i.e., we multiply the contin- 
uum model with 



CYCLABS(£) = exp - 



r(WE/E cyc ) 2 
(E - E cyc ) 2 + W 2 



(7) 



with the centroid energy, E cyc , the width of the feature, W, and 
the optical depth of the line, r. This approach allows the com- 
parison of the line parameters with results from previo us papers, 
in which also the CYCLABS model has been used (e.g., [Nakajima 
et al.|2006[rTsygankov et al.|2007[|Li et al.|2012| ). In the remain- 
der of the paper, we label the cyclotron line parameters with the 
number of respective harmonic, where corresponds to the fun- 
damental line. As common for CRSFs, W and r are strongly cor- 
related ( |Coburn|2001| ). We therefore set a lower limit of 0.5 keV 
for the width, which is comparable to the typical PCA resolution. 

In order to describe the data of 4U 01 15+634, in addition to 
several of the continuum models discussed above, we take in- 
terstellar absorption into account, modelling it with an upd ated 
version ofTBabfl using abundances by | Wilms et al.| ( |2000| ) and 
cross sections by |Verner & Yakovl"ev| ( |1995| ). We furthermore ac- 
count for an intrinsic Fe Ka flourescence line, which we describe 
with an additive narrow Gaussian emission line with frozen cen- 
troid energy at 6.4 keV, a width of 10" 4 keV, and a flux A Fe and 
equivalent width W Fe . Assuming a narrow line is justified be- 
cause in systems like 4U 0115+63 Fe Ka is expected to have a 
width in the order of <0.5 keV (see, e.g., Torrejon et al.|2010| ), 
which cannot be resolved by PCA. 

In order to take into account that the flux normalization of the 
different instruments used is not perfectly known we introduce 
cross calibration constants chexte and cisgri as free fit param- 
eters. These constants also account for flux differences of the 
source between the observations which were not fully simulta- 
neous. The PCA background was estimated using the standard 
models provided by GSFC. These estimates often show a flux 
deviation from the proper background by a few percent. This de- 
viation was accounted for by multiplying the background flux 



with another fit parameter, (see also Rothsc hild et al. 201 1 
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Fig. 2. Deviations between the data and the best fit with 
Cutoff PL and lOkeV feature for the PCA (circles), HEXTE 
(crosses), and ISGRI (filled diamonds) in units of standard devi- 
ations o~. For these fits, no cyclotron lines have been taken into 
account. The numbers on the right side indicate the respective 
epoch. 



3.2. Fitting Strategy 

Since the Cut of fPL is the easiest continuum model for X-ray 
pulsars, we used this model for our fits. Furthermore, the lOkeV 
feature in emission is required to get a good description of the 
data. For each epoch, the fits have been simultaneously per- 
formed for all available intstruments. Figure [2] shows the residu- 
als of the best fit using this continuum without taking any CRSFs 
into account. These residuals clearly indicate the need of the 
cyclotron lines in this model at ~11, 22, and sometimes above 
30 keV. As the quality of individual data sets is strongly variable 
over the outburst, in an initial fit run a variable number of CRSFs 
were added to individual fits until the reduced^ 2 , ^ 2 ed ~ 1. To 
avoid continuum modeling with the absorption features, we fixed 
the widths of the CRSFs of the second and higher harmonics to 
typical values, i.e., 4keV (Ferrig no et al.|2 009). 
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Fig. 3. a, c, and d: Confidence contours for the folding energy, 
the photon index and the hydrogen column density for epoch 6. 
b: Confidence contours between the Efoid and T with Nu frozen 
to 1.7 x 10 22 cm" 2 . Contour lines correspond to the 68.3%, 90%, 
and 99% level. Color indicates Ax 2 with respect to the best fit 
value, with the color scale running from orange (low A^ 2 ) to 
dark blue (large A^ 2 ). 
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Fig. 4. a: Folded spectrum from epoch 6 together with best fit 
(red) and all model components, i.e., Cutoff PL, lOkeV fea- 
ture, Fe Ka flourescence line, and five cyclotron features, b: 
Residuals of the best fit, shown as the difference between the 
data and the model normalized to the uncertainty of each data 
point. Dark blue data points correspond to PCA, cyan points to 
HEXTE, and orange points to ISGRI. 
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Fig. 5. Left Temporal evolution and right luminosity depen- 
dence of the continuum parameters (see text for definitions) and 
lOkeV feature for the fits with the Cut off PL continuum. The 
colour gradient from reddish to blueish data points indicates the 
temporal evolution of the outburst. The data points show results 
from the initial fit run, where a variable number of cyclotron 
lines was used. Green data points in the left column correspond 
to the second fit run with a fixed number of three cyclotron lines 
(only in the left column), which was been performed as a consis- 
tency check (see text for details). Where no green points are vis- 
ible, both fits gave the same results. Lx is in units of 10 37 erg s" 1 . 
Aq and Ap e are in units of photons s" 1 cm -2 , while Ap e is in mul- 
tiples of 10" 3 . 



In order to check whether the number of different cyclotron 
features in the model affects the fit results of the continuum pa- 
rameters and the fundamental and first harmonic CRSF, in a sec- 
ond run all spectral fits included the fundamental line and its first 
and second harmonic only. For epochs, where only the first har- 
monic was necessary to obtain a good fit in our initial run, we 
added the second harmonic with fixed parameters during the sec- 



ond run, holding its parameters at the mean values of the results 
from the other epochs. 

Since initial fits with the Cut of fPL continuum showed large 
uncertainties for the hydrogen column density, TVh, without a sig- 
nificant variability between different epochs, we held A^h fixed at 
the mean value A^h = 1.7 X 10 22 cm -2 to avoid unphysical corre- 
lations between these fit parameters (Fig. [3]). 
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Fig. 6. Luminosity and time dependence (inset) of the energy of 
the fundamental CRSF using the Cutoff PL continuum model. 
The color coding is the same as in Fig. [5] Cyan data points (dia- 
monds) show the results when using the NPEX continuum. 



The two fit runs show that almost all fit parameters are equal 
within their respective uncertainties for both runs. Only £f old and 
the line parameters of the second harmonic are slightly different 
in epoch 6 and 7. These differences can be explained by the fact 
that for these two cases the second run has to account for the 
equivalent of the four harmonics required in the first run. The 
differences affect mainly the higher energies and are compen- 
sated by a change in the folding energy and the closest cyclotron 
line, i.e., the second harmonic. The remainder of this paper is 
based on the results from the first fit run. 

An example for the best fit using the Cut off PL contin- 
uum with lOkeV feature and five cyclotron features (fit run a, 
epoch 6) is shown together with all model components in Fig. 4] 
This spectrum was recorded by all three instruments during the 
maximum of the outburst, i.e., epoch 6. All fit results from runs 
using the Cutoff PL continuum are summarized in Table [2] 

4. Spectral Evolution: No Anti-Correlation of the 
Cyclotron Line Energy with Flux 



Figures [5}j7] show the evolution of the continuum and cyclotron 
line parameters as well as the parameters for the Fe Ka line for 
the the first fit run with the Cutoff PL continuum and lOkeV 
feature, and fixed Nn (see also Tab. [2]). As already noted in pre- 
vious analyses for 4U 0115+634, the continuum varies over the 
outburst (e.g., |Li et al.|2012| ). During the brightest phases of the 
outburst the spectrum is the hardest one. As shown in Fig. [3] the 
variation of T is significantly larger than the very slight corre- 
lation seen between T and £f id- As expected, the source's in- 
trinsic iron Ka line is significantly present during all epochs. As 
also observed in other X-ray transient pulsars ( |Inoue|1985] ) and 
expected for a line originating from fluorescence, the line flux 
is positively correlated to the X-ray flux from the source. The 
equivalent width stays constant within the uncertainties over the 
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outburst. The lOkeV feature varies in both, the centroid energy 
E G and its flux A G . Its energy seems to be slightly anti-correlated 
with the source flux, while the relative flux in the feature, i.e., 
the ratio between the flux in the the 10 keV feature and the total 
flux, remains constant. Its width, cr G , remains relatively constant 
between 3.0 keV and 3.5 keV. 

Very different from virtually all previous analyses, however, 
the cyclotron line is seen to be only slightly varying over the out- 
burst. The colored data points in Fig. [5] and [7] show the centroid 
energy, Eq, the width, Wo, and the optical depth, To, of the fun- 
damental cyclotron line against time and flux for the Cutoff PL 
continuum, modified with the lOkeV feature. In these fits, the 
width Wo remains constant throughout the outburst, and there is 
only a slight (~3cr) indication that the optical depth To was some- 
what shallower during the initial phases of the outburst, before 
~MJD 54552, although it cannot be excluded that some of this 
effect is due to a correlation with the "lOkeV feature", which 
during this time of the outburst peaks close to the cyclotron 
line at around 8keV. Most importantly, however, as shown in 
Fig. [6] the energy of the fundamental cyclotron line is only very 
slightly variable: during early phases of the outburst the line 
is around llkeV, moving towards lOkeV as the outburst pro- 
gressed. This result is contrary to all earlier studies of outbursts 
of 4U 0115+634, where strong changes in the line energy were 
seen (e.g.,|Li et al.|2012|JTsygankov et al.|2007[|Nakajima et al. 



2006; Mi hara et ali|2004] 



Where does this difference in line behavior come from? 



As discussed in Sect. 3.1 a large variety of continuum mod- 
els can be used to describe the broad band spectra of accret- 
ing neutron stars. Due to the early successes of fitting spectra of 
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Table 2. Continuum parameters of the time resolved spectral analysis. Fits were performed with the Cutoff PL model and a varying 
number of CRSFs. 
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Notes. Uncertainties and upper limits are at the 90% confidence level for one interesting parameter. 

{a) Epoch for data grouping (see tex t for details). Absorbed lum inosity in units of 10 37 erg s" 1 . L x covers the energy band 3-50 keV and was 
calculated using a distance of 7 kpc ( |Negueruela & O kazaki 2001). (c) In units of photons s _1 cm -2 . 



4U 0115+634 with the NPEX and PLCUT models, these continua 
have been extensively used when modeling the outburst behavior 
of this pulsar. Unfortunately, however, in most previous analyses 
- including some of our own - no attempt was made at modeling 
the data with any of the other available models. 

As an example, modeling the 2008 outburst of 4U 1 15+634 
with the NPEX model (following Nakajima et al. 2006 , fixing T2 
to 2.0) results in fits that are of only slight worse quality as those 
using the Cutoff PL model modified by the lOkeV feature, and 
we can recover the strong variation of the cyclotron line found 
in earlier analyses (Fig. [6] cyan data points^] Depending on the 
choice of the continuum, we therefore find a fundamentally dif- 
ferent behavior of the fundamental cyclotron line, especially the 
centroid energy. 



4 Modeling the spectra with the PLCUT model results in an unphysical 
break around the energy of the fundamental line, we therefore consider 
this continuum to be unsuitable for studying the cyclotron line behavior. 



Given that the fits with NPEX and Cut of fPL give similarly 
good results, what is the fundamental difference between both 
modeling approaches? NPEX models generally result in very 
broad cyclotron lines, with Wp and W\ often exceeding 5keV 
(e.g., [Nakajima et al. 2006 runs using the PLCUT model give 
similar results, see, e.g., |Li et al.||2012| . In contrast, using the 
Cutoff PL continuum with the lOkeV feature results in widths 
of typically less than 3 keV, much more consistent with the nar- 
row shapes of the residuals shown in Fig. [2] When lines are as 
broad as in the NPEX fits, they can influence the continuum fit. 
To illustrate how strongly these features actually distort the con- 
tinuum, Fig. [8] shows the continuum shape inferred when setting 
To and ri to zero. The residuals for the CutoffPL continuum 
with the 10 keV feature (dark blue) show sharp absorption dips at 
the cyclotron line energies of the fundamental and first harmonic 
line energies. For the NPEX model (green line), on the other hand, 
it is obvious that rather than describing narrow cyclotron lines, 
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cyclotron lines to zero, b Ratio between the data and the model. 



the multiplicative broad line model strongly influences the con- 
tinuum. As discussed above, our spectral fits show that during 
the brightest phases of the outburst a significant hardening of the 
continuum emission is observed, together with a change of the 
exponential cutoff. While this is also seen in the variation of the 
continuum parameters, it is very likely that the luminosity de- 
pendence of the CRSF in the NPEX fits is due to the line partially 
modeling this behavior of the continuum, and not a real physical 
effect. 

Using the results from the CutoffPL plus lOkeV feature fits 
for the line energy also solves another problem found in ear- 
lier spectral modeling. Here, using NPEX fits, the ratios between 
the fundamental and higher cyclotron lines were often found to 
deviate significantly from integers. While slightly non-integer 
ratios are expected when taking relativistic quantum mechan- 
ics into account (see the discussion by Pottschn iidt et al.|2005| ), 
the large deviations from integer multiples seen previously (e.g., 



Tsygankov etal.|2007; 


Nakajima et al.|2006 JHeindl et al.|1999| 


Santangelo et al.||1999 


) require rather complex model assump- 



tions which introduce second order effects such as strong vertical 
5-field gradients, crustal field structures on small scales, or ther- 



momagnetic effects for their explanation (see Schon herr et al. 
|2007| and references therein for a discussion). Figure [9] shows 
that we can recover these non-integer ratios in our NPEX fits, but 
when using the CutoffPL plus lOkeV feature model, these ra- 
tios are mostly in agreement with integer values, or at least only 
slightly higher, as expected from the most simple models for cy- 
clotron line formation. We conclude that also from a physical 
point of view, the CutoffPL continuum with the lOkeV feature 
gives a much more satisfactory description of the data. 

Finally, we check the dependence of the line position on 
other free fit parameters when using the CutoffPL continuum 
modified by the lOkeV feature. Since the origin of the lOkeV 
feature remains still an open questtion, possible relations be- 
tween the CRSF centroid energy and the parameters of the 
lOkeV feature have to be investigated, and other continuum pa- 




60 70 
time (MJD-54500) 

Fig. 9. Ratios of the cyclotron line energies E n /Eo with respect 
to the centroid energy of the fundamental line against time. The 
dark green, dark blue, red and purple data points (bars) corre- 
spond to E\/Eo, E2/E0, E3/E0, and E4/E0, respectively, using 
the CutoffPL continuum with a lOkeV feature. The light green 
and cyan data points (filled circles) show the results for E\/Eq 
and E2/E0 using NPEX. The dashed lines indicate integer values 
for these ratios. 



rameters such as Ef \ d or T could influence the results for the cy- 



lotron line parameters. Figure 10 shows the behavior of A% in 



the vicinity of the best fit values for three epochs representative 
for the early phase, peak, and late phase of the 2008 outburst for 
selected combinations of important fit parameters. These con- 
tour plots show that the line energy is not affected by cross cor- 
relations with other parameters, i.e., possible variations in the 
line energy of the fundamental cyclotron line are not due to vari- 
ations observed in the parameters of the lOkeV feature or the 
continuum parameters. 

In contrast to the stable behavior of the line when using 
CutoffPL with a lOkeV feature, Fig. \TT\ illustrates the depen- 
dence of the CRSF centroid energy on the continuum parameters 
of the NPEX model. In many cases, strong correlations between 
Eq and T\ as well as between Eq and £f id are present. These cor- 
relations are a further indicator that the Lorentzian component is 
used to model the continuum rather than a cyclotron line when 
using the NPEX model. 



5. Summary and Conclusion 

In this paper we have presented a study of the 2008 outburst 
of 4U 01 15+634 based on RXTE and INTEGRAL data (Fig. (I). 
We reproduced the results from previous work (e.g., [Nakajima 
|eTaL|2006l |Tsygankov eT^2WJ\ |Li et al.||20T2l and others), 
that the spectra can be successfully modelled by, e.g., the NPEX 
model. We showed that for these fits the very broad absorption 
features, which are thought to describe the CRSFs, rather model 
the broadband continuum. Our result that the continuum param- 
eters are strongly variable over the outburst (Fig. [5]), could be 
responsible for the change in the cyclotron line energy in these 
fits. 



8 



Sebastian Miiller et al.: No anti-correlation between cyclotron line energy and X-ray flux in 4U 01 15+634 




10 12 14 

#o [keV] 



10 12 14 

#o [keV] 



10 12 14 

#o [keV] 



Fig. 10. Confidence contours between the fundamental CRSF 
energy and other fit parameters. The columns show the results 
for epochs 2, 6, and 14, respectively. Contour lines correspond 
to the 68.3%, 90%, and 99% confidence level. Color indicates 
Ax 2 with respect to the best fit value, with the color scale run- 
ning from orange (low Ax 2 ) to dark blue (large Ax 2 ). 



We have shown that the spectral continuu m can be also well 
described with a model introduced by Klochkov et al.| ( [2007] ) for 
EXO 2030+375 and |Ferrigno eTal] ( |2009) for 4U 0115+634, 
i.e., the Cutoff PL, modified by strong Gaussian emission fea- 
ture around lOkeV and several (up to four) cyclotron lines 
(Tableland Fig.|4j. Due to an unphysical break close to the en- 
ergy of the fundamental cyclotron line, the PLCUT model does 
not give a physical description of the continuum. Using the 
Cutoff PL together with the lOkeV feature, we have shown that 
the modeling of lines reflects the theoretical expectation in both 
the line shape and the ratios of the harmonics. In this model, the 
luminosity dependency of the centroid energy is not present. As 
shown also by Ferri gno et al.| ( |20lT] ), the line centroid energy is 
remarkably stable throughout the brightest phase of the outburst. 



15 
^ 10 

b 5 

3 
2 



> 



+ 


L + 






aMiMK. j jam 










+ 







10 12 14 

E [keV] 



10 12 14 

E [keV] 



10 12 14 

E [keV] 



Fig. 11. Confidence contours between the fundamental cyclotron 
line energy and continuum parameters when using the NPEX 
model. The displayed epochs, colours and lines are the same as 
in Fig.fTOl 



Here we show that such stability is present, albeit at lower signif- 
icance, also during the early and late phases. In this model, the 
continuum variations are explained by a slight anti-correlation 
of the centroid energy of the Gaussian feature with luminosity 
and a significant variation of T. 

The strong systematic influence of the chosen continuum on 
the cyclotron line behavior illustrates a problem that is poten- 
tially inherent to all cyclotron line measurements. We empha- 
size, however, that this does not mean that the general picture 
of the different regimes of magnetized neutron star accretion 
and gene ral luminosity depen dence of the cyclotron line as out- 
lined by Becke Fet al.| ( |2012| and references therein) is incor- 
rect. In the sample of sources discussed by Bec ker et al. ( 2012 ), 
4U 0115+634 was an outlier. The change of line energy inferred 
from the NPEX and PLCUT models was almost a jump at a 3- 
50keV luminosity of around 3 x 10 37 ergs _1 and strong hys- 
teresis effects were present (Fig. [6]). Such a behavior has not 
been seen in any of the other cyclotron line sources. The sec- 
ond strongest luminosity-depen dent cyclotron line variability is 
observed in V0332+53 (Mowlavi et al.||2006[ |Tsygankov et al" 



2006; Nakaji ma et al.|2010| ), in which the cyclotron lines are so 
strong that the change can be seen by eye even in the raw detec- 
tor spectra. While the boundaries between the different accretion 
regimes estimated by |Becker et al.| ( |2012| ) might therefore have 
to be adjusted, the overall principal behavior of the accretion 
column discussed by these authors still holds. 



With the recently successful launch of NuSTAR (Harrison 
|et al.|2010| ), a low background, focusing telescope with slightly 
higher resolution than the RXTE-HEXTE or the Suzaku-HXD is 
now available. With its higher resolution and the much higher 
expected signal to noise ratio in the spectra, we will be able, 
for the first time, to resolve cyclotron lines and study their 
luminosity-dependent behavior without the danger of a system- 
atic influence of the choice of a continuum. 



9 



Sebastian Miiller et al.: No anti-correlation between cyclotron line energy and X-ray flux in 4U 01 15+634 



Acknowledgements. We thank the referee for his/her insightful comments. 
We also thank the schedulers of RXTE and INTEGRAL for their role in 
making this campaign possible, and the International Space Science Institute 
in Bern, Switzerland, for their hospitality. We acknowledge funding by the 
Bundesministerium fiir Wirtschaft und Technologie under Deutsches Zentrum 
fur Luft- und Raumfahrt grants 50OR0808, 50OR0905, and 50OR1113, and by 
the Deutscher Akademischer Austauschdienst. MTW is supported by the US 
Office of Naval Research. IC acknowledges financial support from the French 
Space Agency CNES through CNRS. SMN and JMT acknowledge support from 
the Spanish Ministerio de Ciencia, Tecnologia e Innovacion (MCINN) through 
grant AYA201 0-1 5431 and the use of the computer facilities made available 
through the grant AIB2010DE-00057. This research is in part based on observa- 
tions with INTEGRAL, an ESA project with instruments and science data centre 
funded by ESA member states (especially the PI countries: Denmark, France, 
Germany, Italy, Switzerland, Spain), Czech Republic, and Poland, and with the 
participation of Russia and the USA. We thank John E. Davis for the develop- 
ment of the SLxf ig module, which was used to create all figures in the paper. 



References 

Arnaud, K. A. 1996, in Astronomical Data Analysis Software and Systems V, 
ed. J. H. Jacoby & J. Barnes, ASP Conf. Ser. 101 (San Francisco: Astron. 
Soc. Pacific), 17 
Basko, M. M. & Sunyaev, R. A. 1976, MNRAS, 175, 395 
Becker, P. A., Klochkov, D., Schonherr, G., et al. 2012, A&A, 544, A123 
Becker, P. A. & Wolff, M. T. 2007, ApJ, 654, 435 
Bradt, H. V., Rothschild, R. E., & Swank, J. H. 1993, A&AS, 97, 355 
Burnard, D. J., Arons, J., & Klein, R. I. 1991, ApJ, 367, 575 
Caballero, L, Kretschmar, P., Santangelo, A., et al. 2007, A&A, 465, L21 
Caballero, I. & Wilms, J. 2011, Mem. Soc. Astron. Ital., 83, 230 
Coburn, W. 2001, Ph.D. Thesis, University of California, San Diego 
Cominsky, L., Clark, G. W, Li, E, Mayer, W., & Rappaport, S. 1978, Nat, 273, 
367 

Davidson, K. 1973, Nat, 246, 1 

Ferrigno, C, Becker, P. A., Segreto, A., Mineo, T., & Santangelo, A. 2009, A&A, 
498, 825 

Ferrigno, C, Falanga, M., Bozzo, E., et al. 2011, A&A, 532, A76 

Giacconi, R., Murray, S., Gursky, H., et al. 1972, ApJ, 178, 281 

Harrison, F. A., Boggs, S., Christensen, E, et al. 2010, in Space Telescopes 

and Instrumentation 2010: Ultraviolet to Gamma Ray, ed. M. Arnaud, S. S. 

Murray, & T. Takahashi, SPIE Conf. Ser. 7732, 77320S 
Heindl, W. A., Coburn, W., Gruber, D. E., et al. 1999, ApJ, 521, L49 
Houck, J. C. & Denicola, L. A. 2000, in Astronomical Data Analysis Software 

and Systems IX, ed. N. Manset, C. Veillet, & D. Crabtree, ASP Conf. Ser. 216, 

591 

Inoue, H. 1985, Space Sci. Rev, 40, 317 

Jahoda, K., Markwardt, C. B., Radeva, Y., et al. 2006, ApJS, 163, 401 
Johns, M., Koski, A., Canizares, C, et al. 1978, IAU Circ. 3171 
Klochkov, D., Horns, D., Santangelo, A., et al. 2007, A&A, 464, L45 
Klochkov, D., Santangelo, A., Staubert, R., & Ferrigno, C. 2008, A&A, 491, 833 
Klochkov, D., Staubert, R., Santangelo, A., Rothschild, R. E., & Ferrigno, C. 

2011, A&A, 532, A126 
Kretschmar, P., Kreykenbohm, I., Wilms, J., et al. 1997, in The Transparent 

Universe, ed. C. Winkler, T. J.-L. Courvoisier, & P. Durouchoux, ESA Special 

Publication 382, 141 
Kreykenbohm, I., Coburn, W., Wilms, J., et al. 2002, A&A, 395, 129 
Kreykenbohm, I., Kretschmar, P., Wilms, J., et al. 1999, A&A, 341, 141 
Krimm, H. A., Barthelmy, S. D., Baumgartner, W, et al. 2008, Astron. Tel. 1427 
Lebrun, E, Leray, J. P., Lavocat, P., et al. 2003, A&A, 411, L141 
Li, J., Wang, W, & Zhao, Y. 2012, MNRAS, 423, 2854 
Makishima, K., Mihara, T., Nagase, E, & Tanaka, Y. 1999, ApJ, 525, 978 
Meszaros, P. 1984, Space Sci. Rev., 38, 325 

Meszaros, P. 1992, High-energy radiation from magnetized neutron stars 

(Chicago: Chicago Univ. Press) 
Mihara, T., Makishima, K., & Nagase, F. 2004, ApJ, 610, 390 
Mihara, T., Makishima, K., Ohashi, T., Sakao, T., & Tashiro, M. 1990, Nat, 346, 

250 

Mowlavi, N., Kreykenbohm, I., Shaw, S., et al. 2006, A&A, 451, 817 

Miiller, S., Obst, M., Kreykenbohm, I., et al. 2010, in proceedings of "Eighth 

Integral Workshop. The Restless Gamma-ray Universe (INTEGRAL 2010)", 

PoS(INTEGRAL 2010)1 16 
Nagase, E, Dotani, T., Tanaka, Y, et al. 1991, ApJ, 375, L49 
Nakajima, M., Mihara, T., & Makishima, K. 2010, ApJ, 710, 1755 
Nakajima, M., Mihara, T., Makishima, K., & Niko, H. 2006, ApJ, 646, 1125 
Negueruela, I. & Okazaki, A. T. 2001, A&A, 369, 108 
Pottschmidt, K., Kreykenbohm, I., Wilms, J., et al. 2005, ApJ, 634, L97 



Rappaport, S., Clark, G. W, Cominsky, L., Li, E, & Joss, P. C. 1978, ApJ, 224, 
LI 

Rothschild, R. E., Blanco, P. R., Gruber, D. E., et al. 1998, ApJ, 496, 538 
Rothschild, R. E., Markowitz, A., Rivers, E., et al. 2011, ApJ, 733, 23 
Santangelo, A., Segreto, A., Giarrusso, S., et al. 1999, ApJ, 523, L85 
Schonherr, G., Wilms, J., Kretschmar, P., et al. 2007, A&A, 472, 353 
Schwarm, F. W, Becker, P. A., Schonherr, G., et al. 2012, A&A, in prep. 
Staubert, R., Shakura, N. I., Postnov, K., et al. 2007, A&A, 465, L25 
Tanaka, Y. 1986, in Radiation Hydrodynamics in Stars and Compact Objects, ed. 

D. Mihalas & K. H. Winkler, IAU Coll. 86, New York, Springer, 198 
Torrejon, J. M., Schulz, N. S., Nowak, M. A., & Kallman, T. R. 2010, ApJ, 715, 

947 

Tsygankov, S. S., Lutovinov, A. A., Churazov, E. M., & Sunyaev, R. A. 2006, 
MNRAS, 371,19 

Tsygankov, S. S., Lutovinov, A. A., Churazov, E. M., & Sunyaev, R. A. 2007, 

Astronomy Letters, 33, 368 
Unger, S. J., Roche, P., Negueruela, I., et al. 1998, A&A, 336, 960 
Verner, D. A. & Yakovlev, D. G. 1995, A&AS, 109, 125 
Wheaton, W. A., Doty, J. P., Primini, F. A., et al. 1979, Nat, 282, 240 
White, N. E., Swank, J. H., & Holt, S. S. 1983, ApJ, 270, 711 
Wilms, J., Allen, A., & McCray, R. 2000, ApJ, 542, 914 
Winkler, C, Courvoisier, T. J.-L., Di Cocco, G., et al. 2003, A&A, 411, LI 



10 



